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O ■ Abstract 
O. 

The dynamically generated effective gluon mass is known to depend non-trivially on the mo 
mentum, decreasing sufficiently fast in the deep ultraviolet, in order for the renormalizability of 
QCD to be preserved. General arguments based on the analogy with the constituent quark masses, 
as well as explicit calculations using the operator-product expansion, suggest that the gluon mass 
falls off as the inverse square of the momentum, relating it to the gauge-invariant gluon condensate 
of dimension four. In this article we demonstrate that the power-law running of the effective gluon 



Oh 



X 



mass is indeed dynamically realized at the level of the non-perturbative Schwinger-Dyson equation. 
We study a gauge-invariant non-linear integral equation involving the gluon self-energy, and estab- 



lish the conditions necessary for the existence of infrared finite solutions, described in terms of a 
momentum-dependent gluon mass. Assuming a simplified form for the gluon propagator, we derive 
a secondary integral equation that controls the running of the mass in the deep ultraviolet. De- 

00 i 

pending on the values chosen for certain parameters entering into the Ansatz for the fully-dressed 
three-gluon vertex, this latter equation yields either logarithmic solutions, familiar from previous 



linear studies, or a new type of solutions, displaying power-law running. In addition, it furnishes 
a non-trivial integral constraint, which restricts significantly (but does not determine fully) the 
running of the mass in the intermediate and infrared regimes. The numerical analysis presented is 
in complete agreement with the analytic results obtained, showing clearly the appearance of the 
two types of momentum-dependence, well-separated in the relevant space of parameters. Several 
technical improvements, various open issues, and possible future directions, are briefly discussed. 

PACS numbers: 12.38.Lg, 12.38.Aw 
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I. INTRODUCTION 



The possibility that the non-perturbative dynamics of QCD generate an effective gluon 

, has received considerab 



mass, first elaborated in the pioneering work of Cornwal 
theoretical g fl flfl, 0, ^dphe—logicalF 



1CJ, attention in recent years [12]. 



As advocated in 2|, the origin of this effective mass is purely dynamical and preserves 
the local SU(3) C invariance of QCD, in close analogy to what happens in QED 2 (Schwinger 
model) [13I , where the photon acquires a mass without violating the abelian gauge symmetry. 
The gluon mass is not a directly measurable quantity, and its value is determined, at least 
in principle, by relating it to other dimensionful non-perturbative parameters, such as the 



string tension, glueball masses, gluon condensates, and the vacuum energy of QCD 



Since gluon mass generation is a purely non-perturbative effect, the most standard way for 



studying it in the continuum is through the Schwinger- Dyson equations (SDE) governing t 
relevant Green's functions. These equations have been treated from the very beginning 



lie 

21 



in a manifestly gauge-invariant way, by resorting to 



re pinch technique (PT) 
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;he systematic rearrangement of graphs 
Q. Subsequently, the powerful all- 



implemented by t 

order connection 16fl between the PT and the Feynman gauge of the Background Field 

rT nn 

Method (BFM) [17| gave rise to the "PT-BFM" truncation scheme [8|, [l8| that guarantees 
crucial properties, such as gauge-invariance, gauge-independence, and invariance under the 
renormalization-group (RG). What one characterizes as gluon mass generation at the level 
of the SDE governing the PT-BFM gluon self-energy is essentially the existence of solutions 
that reach a finite (non- vanishing) value in the deep infrared (IR). These solutions may be 
successfully fitted by a "massive" propagator of the form A _1 (g 2 ) = q 2 + m 2 (q 2 ); the crucial 
characteristic, enforced by the SDE itself, is that m 2 (q 2 ) is not "hard", but depends non- 
trivially on the momentum transfer q 2 . Specifically, m 2 (q 2 ) is a monotonically decreasing 
function, starting at a non- zero value in the IR (m 2 (0) > 0) and dropping "sufficiently fast" 
in the deep ultraviolet (UV). When the RG logarithms are properly taken into account, one 
obtains in addition the non-perturbative generalization of g 2 {q 2 ), the QCD running coupling 
(effective charge). The presence of m 2 (q 2 ) in the corresponding logarithms tames the Landau 
singularity associated with the perturbative f3 function, and the resul ting effective c harge i s 
asymptotically free in the UV, "freezing" at a finite value in the IR 
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23]. 



The running of m 2 (q 2 ) is of central importance for the self-consistency of this entire 
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approach. Roughly speaking, the value of A _1 (0) is determined by integrals involving A(g 2 ), 
m 2 (q 2 ), and g 2 (q 2 ) over the entire range of (Euclidean) momenta. The UV convergence of 
these integrals depends crucially on how m 2 (q 2 ) behaves as q 2 — > oo. If m 2 (q 2 ) drops off 
asymptotically faster than a logarithm, then A _1 (0) is finite. This is crucial because the 
finiteness of A _1 (0) guarantees essentially the renormalizability of QCD. Had the mass been 
constant instead, these integrals would diverge quadratically; to absorb such a divergence 
one would have to introduce a counterterm of the form ■rriQ(A 2 fV )A 2 at the level of the 
fundamental QCD Lagrangian, which is clearly forbidden by the local gauge invariance. 

The UV behavior of m 2 (q 2 ) is not imposed by hand, but is instead determined dynam- 
ically from the corresponding SDE. The situation is conceptually very similar to the case 
of the dynamically generated (constituent) quark masses, whose momentum dependence is 
controlled by the gap equation. Notice, however, that there is an important technical differ- 
ence. Due to its Dirac structure the gap equation may be separated into two independent 
components, one determining the wave function and one the mass of the quark self-energy. 
In the case of the SDE for the gluon self-energy there is no such direct separation: instead 
an appropriate matching of the contributions on both sides of the equations must be carried 
out. 

In previous studies of linear SDE equations 0, Q] the UV running of m 2 (q 2 ) has been 
found to be logarithmic, of the general form m 2 (q 2 ) ~ (lng 2 ) -1 " 7 , with 7 > in order 
for the aforementioned integrals determining A _1 (0) to converge. Of course, it is natural 
to ask whether the QCD dynamics allows also for a m 2 (q 2 ) displaying power-law running, 



i.e. m 

masses 



q ) ~ q (lng ) 7 , as happens in the case of the dynamically generated quark 



241 ]. This possibility was first envisaged by Cornwall [2j, [25|, based on the afore- 
mentioned studies of chiral symmetry breaking. A decade later, Lavelle wrote down the 
operator-product expansion (OPE) for the (partially dressed) one-loop PT gluon self-energy, 
expressing it solely in terms of the gauge- invariant gluon condensate (G 2 ) = (0| :G°L, G% v : |0) 
of dimension four (no quarks were considered) [23]. The resulting self-energy was then iden- 
tified with an effective gluon mass, m 2 (q 2 ) ~ (G 2 )/q 2 , i.e. a mass displaying power-law 
running. However, to date, the power-law running of the effective gluon mass has never 
been demonstrated as an explicit dynamical possibility at the level of a (non-linear) SDE. 

In this article we show that the non-linear SDE for the gluon self-energy in the PT-BFM 
formalism has two distinct types of solutions: (i) solutions of the type already encountered in 
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the linear studies, with m 2 (q 2 ) running as an inverse power of a logarithm, and (ii) solutions 
found for the first time, where the effective gluon mass drops asymptotically as an inverse 
power of the momentum (multiplied by powers of logarithms). Which of the two types of 
solutions will be actually realized is a complicated dynamical problem, depending mainly 
on the details of the (fully dressed) three-gluon vertex, T^ a/3 , entering into the SDE for the 
gluon self-energy (see Fig. [[]). 

The article is organized as follows: In Sec. [Til after setting up the notation, we present 
the derivation of the SDE for the PT-BFM gluon propagator. We restrict our analysis to 
the subset of "one-loop dressed" gluonic contributions, namely the two fully-dressed gluonic 
diagrams shown in Fig. [TJ In addition, we explain briefly how the naive all-order Ward 
identity (WI) satisfied by the full three-gluon vertex in this formalism enforces the transver- 
sality of the gluon self-energy even in the absence of ghost loops. In Sec. II I II we manipulate 
the SDE derived in the previous section further, with the motivation to search for infrared- 
finite solutions. An important ingredient at this stage is the Ansatz introduced for the full 
three gluon vertex, based on the gauge-technique; this vertex satisfies, by construction, the 
simple all-order WI characteristic of the PT-BFM, and contains sufficient structure to give 
rise to a non- vanishing A _1 (0). In addition, we discuss in detail the technical adjustments 
implemented to the SDE, in order to endow it with the correct RG behavior. The next two 
sections contain the main results of this article. Specifically, in Sec. IIVI we extract from the 
SDE an integral equation that determines the running of the effective gluon mass in the UV, 
and study its solutions. We derive an important constraint, in the form of an integral bound- 
ary condition, relating the value of A _1 (0) with m 2 (q 2 ) in the entire range of momenta; this 
condition must be necessarily satisfied in order for the mass equation to have solutions that 
vanish asymptotically in the deep UV. Then we demonstrate that, depending on the values 
of two basic parameters, one finds solutions for the masses that drop as inverse powers of a 
logarithm of q 2 , or much faster, as an inverse power of q 2 . In Sec. |V]we carry out a numerical 
analysis of the SDE, which fully corroborates the existence of the two aforementioned types 
of solutions. In Sec. IVII we present our conclusions and a discussion of various open issues. 
Finally, in an Appendix we present all technical points related to the modifications one must 
introduce to the standard angular approximation in order to capture correctly the leading 
m 2 (q 2 ) behavior, encoded in the exact SDE. 
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II. SDE FOR THE PT-BFM GLUON PROPAGATOR 



In this section we derive in detail a non-linear SDE equation for the gluon propagator in 
the PT-BFM formalism. As has been explained in the literature, this scheme is essentially 
founded on the all-order correspondence between PT and BFM [3]: the (gauge- independent) 
PT effective n-point functions coincide with the (gauge-dependent) BFM n-point functions 
provided that the latter are computed in the Feynman gauge. One of the most powerful 
features of this formalism is the special way in which the transversality of the all-order 
PT-BFM self-energy is realized. Specifically, by virtue of the Abelian-like Wis satisfied by 
the vertices involved, gluonic and ghost contributions are separately transverse, within each 
order in the "dressed-loop" expansion of the corresponding SDE [8]. This property, in turn, 
allows for a systematic truncation of the full SDE 



that preserves the crucial property 



of gauge invariance. In particular, instead of a system of two coupled equations involving 
gluon and ghost propagators, one may consider only the subset containing gluons, without 
compromising the transversality of the gluon self-energy. Therefore, in what follows, we will 
consider only the gauge-invariant subset of "one-loop dressed" gluonic diagrams, given by 
the two graphs of Fig. [TJ 

Within this formalism there are two distinct gluon propagators, the background gluon 
propagator A fJil/ (q) ( which, in the Feynman gauge, coincides with the PT gluon propaga- 
tor) and the quantum gluo n p ropagator A^ u (q), appearing inside the loops. By virtue of a 



powerful all-order identity 29. 130 1 . one may express A Mi ,(g) in terms of A^ u (q) and auxiliary 



(unphysical) Green's functions involving anti-fields and background sources. As a first ap- 
proximation, in this work we will neglect the effects of the aforementioned auxiliary Green's 
functions, and carry out the substitution A fMU (q) — > A fll/ (q) throughout (and drop the "hats" 
to simplify the notation). 

In the Feynman gauge (both of the usual linear renormalizable gauges as well as the 
BFM) the full gluon propagator A^ v (q) has the general form (we suppress color indices) 



where 



P^(<?)A(g 2 ) + ^ 



(2.1) 



PM= 9^- qj f, (2-2) 
is the usual transverse projector. The scalar function A(g 2 ) is related to the all-order gluon 
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self-energy U^q), 

U llu (q) = P flu (q)U(q 2 ), (2.3) 

through 

A-V) = q 2 + iTl(q 2 ). (2.4) 

Notice that, since n(g 2 ) has been defined in (12. 4p with the imaginary factor i pulled out 

in front, it is given simply by the corresponding Feynman diagrams in Minkowski space. 
Finally, the inverse of the full gluon propagator is given by 

A ^(<?) = WnuQ 2 ~ n^(g) 

= iP M „( g )A-V)+ ( 2 - 5 ) 



At tree- level, 

a™ = 

q 



A^(q) = 9 —, A (q 2 ) = -. (2.6) 



The SDE in Fig. [T] reads, 

(A" 1 )?, (?) = n 2 9^ ah - [n^(?)| ai + K\ a \ > ( 2 - 7 ) 

with the closed expressions corresponding to the diagrams (a\) and (02) given by 
K(<l)\ ai = \ J[dk}f^(q,k,-k-q)A a J(k)^ 

KL = \ fi dk ^£lAZ(k) . (2.8) 

The flow of momenta, together with the color and Lorentz indices, are shown in 
Fig. [TJ We have assumed dimensional regularization, employing the short-hand notation 
[dk] = d d k/ (2ii) d , where d = 4 — e is the dimension of space-time. 

We will use greek letters for the Lorentz indices and latin for the color indices. The 
tree-level vertex appearing in (ai) is given by (pi = k, p 2 = —k — q) 

r^(g,Pi,p 2 ) = gf acx T f , aP (q,p 1 ,p 2 ), 

f m p{q,P\: P2) = (pi - P2)f,g a f5 + 2q p g m - iq a gtf, (2.9) 

and satisfies the elementary WI 

?% a/3 (g,pi,p 2 ) = {p\ -p\)g aP = i [A-^(pi) - A~^(p 2 )] = [Aq 1 ^) - Aq 1 ^)] g a(3 . 

(2.10) 




(ai) 



FIG. 1: The gluonic "one-loop dressed" contributions to the SDE. 



The bare four-gluon vertex T ab £ d appearing in (02) is given by 



r 



abed 
fj,vpa 



f acx f xbd {g^gpa - g^gpv + g w g va ) + f adx f xcb (g pp g va - g pu g pa - g^g pv ) 



+f abx f xcd {g w9va - g pa g pv ) 



■ abc 



(2.11) 



5 ab A„ v and 



In addition, for the fully dressed quantities A ab y and IT ^ we will set A ab u — 

~ abc , — ~ 

= 5 f T ' pa p. In the PT-BFM formalism the all-order WI satisfied by TT^ap is the 
naive generalization of the tree-level WI (12.101) 0, i.e. 

q^^(q, Pl ,p 2 ) = i[A^(pO - A^Cft,)] . (2.12) 

It is then elementary to verify that 



0. 



(2.13) 



Thus, the subset of graphs considered is transverse by itself, i.e. without the inclusion of 
ghost loops, as announced. 

In order to reduce the algebraic complexity of the problem, we drop the longitudinal 
terms from the gluon propagators inside the integrals on the r.h.s. of (12.81) . i.e. we set 
A Q/ g — > —ig a pA. This may be done without compromising the transversality of the answer, 
provided that one drops, at the same time, the longitudinal pieces on the r.h.s. of the WI 
ofEq.flSTH [31|. 

It is then straightforward to arrive at the following form for the SDE 

OV(g)A- V) = ^ ~ ^ (/ l dk \ rfA(k)T uaf3 A{k + q) - 8 g pu j [dk] A(k) 

(2.14) 



where we have used that f ace f bce = S ab CA, with Ca the Casimir eigenvalue in the adjoint 
representation [Ca = N for SU(N)]. After the omission of the longitudinal parts the WI of 
( EZBD becomes 

q v T va p{q, PuP2 )= [A-\p 2 )-A- l ( Pl )]g a p; (2.15) 

contracting by q u and using (12.151) one may easily verify that the sum of the two integrals 
on the r.h.s. of (12. 14j) is indeed transverse. 

III. SDE WITH IR-FINITE SOLUTIONS 

In order to proceed further with Eq. (l2.14p . one needs to supply some information about 
the form of the full vertex IT. To accomplish this, ideally one should set up the corresponding 
SDE governing the vertex T, and solve a system of coupled integral equations. In practice 



this is very difficult and one almost always resorts to the "gauge technique" 32] , expressing 



F as a functional of A, in such a way as to satisfy (by construction) the appropriate WI. 
It is clear that this procedure leaves the transverse (i.e. identically conserved) part of the 
vertex undetermined, a fact that leads to the mishandling of overlapping divergences, and 
forces one to renormalize the resulting SD equation subtractively instead of multiplicatively. 
The Ansatz we will use for the vertex is 



f mP = + Tf* + , (3.1) 



with 



L^(q, Pl ,p 2 ) = r^(q, Pl ,p 2 )+ig a ^ [U( P2 ) - nfa)] 

Zf'Wi.Pa) = -i%(q P g^-q a g^)[Jl{p x )+Jl{p 2 )] , 
q z 



T^{q, PuP2 ) = -ic 2 (qPg» a - q a g^) 



n(pQ n(pa) 

2 ' 2 

Pi PI 



(3.2) 



— fia/3 

Several comments on the properties and role of the vertex IT and its individual com- 



ponents are now in order: 
(i 

(El 



i) When H(pi) = 0, IT^ ^ goes over to the tree-level (bare) result, namely the T^ al3 of 



S 



(ii) 13 is Bose symmetric only with respect to its two legs appearing inside the loop, 
carrying momentum p 1 = k and P2 = — (k + q); so, F M {q,Pi,P2) is invariant under the 
simultaneous exchange pi < — > P2, a < — > [3 and b < — > c. 

(iii) F satisfies by construction the WI of Eq.f l2.12l) . Specifically, when the "longitu- 
dinal" part is contracted with furnishes the r.h.s. of (12.121) . whereas the "transverse" 
parts T± al3 and are identically conserved. 



(iv) F M 13 contains massless poles, i.e. terms going as 1/q 2 ; as is well-known 33|, 



the presence of such poles allows for the possibility A _1 (0) ^ 0. Note that not only the 
longitudinal but also one of the transverse parts contains such poles; this particular feature 
is motivated by the the one-loop analysis of the conventional thrcc-gluon vertex presented 



in 



34|. 



(v) Clearly, the Lorentz structures appearing in (13.11) do not exhaust all tensorial possi- 



bilities. Indeed, the most general parametrization of a vertex with three 
two independent momenta contains 14 linearly independent tensors 



34 



;orentz indices and 



35|; their number 



may be reduced to some extent by imposing Bose symmetry (only partial in our case) and 
the constraint of the corresponding WI. The simplified vertex proposed in (13.11) is only meant 
to capture some of the salient features expected from the full answer, most importantly the 
correct WI, the Bose symmetry, and the presence of poles terms. 

(vi) The transverse parts of F^ ^ are multiplied by the constants C\ and C2, which, at this 
level of approximation, will be treated as arbitrary adjustable parameters. These constants 
affect not only the integral equation but also the two crucial boundary conditions [ (14. 7ft and 
(14.28!) ]. A s we wm see m detail in the next sections, the existence or not of self-consistent 
solutions for the SDE (i.e. solutions satisfying all necessary constraints) depends crucially 
on the values chosen for these constants. At first sight this seems to run contrary to the 
standard lore of the gauge technique, according to which the transverse parts in the Ansatze 
for the vertices do not affect appreciably the solutions in the IR. Notice, however, that, in 
the case we consider here, the distinction between what is of UV and what of IR origin is 
not so sharp. Indeed, as we will see shortly, the value of A _1 (0) - clearly an IR quantity - is 
determined by integrals of the full propagator together with the running mass and coupling 
over the entire range of momenta; there the UV behavior of these quantities is essential. For 
example, the running of the masses in the UV is controlled by the anomalous dimensions, 
which depend themselves on c\ and/or C2- To be sure, in a more complete treatment, the 



values of c\ and c 2 should be uniquely determined by the QCD dynamics, in particular by 
the SDE governing the vertex IT ; unfortunately, such an analysis is beyond our powers 
at the moment. 

We next write (12.141) in the Euclidean space; to that end, we set — q 2 = q 2 , define 
Ae^e) = — A(— q 2 ), and for the integration measure we have [dk] = i[dk] E = id A k E /(27r) 4: . 
Then, using for T ua p the Ansatz of (13.11) . and suppressing the subscript "E", (I2.14p becomes 



A-\q 2 



d 4 k [q 2 + - 



k 2 - 



(k-qf 
q 2 



j A{k)A{k + q) - I d 4 kA(k) 



6bg z 

57T 2 

Qbg' 



J d A k k 2 A(k)A(k + q) - J d 4 kA(k) 
c 2 q 2 J d A k [A(k) - A (k)] A(k + q) , 



(3.3) 



where b = IOCa/^Stt 2 is the contribution of the gluons to the one-loop (3 function in the PT- 
BFM scheme; jthe discrepancy from the correct value b = 11Ca/^8tt 2 is due to the omission 



of the ghosts [8|, [l7| 

The measure in spherical coordinates is given by 

rn poo 

d 4 k = 2n / dx sin 2 \ I dyy ; (3.4) 
Jo Jo 

introducing q 2 = x, k 2 = y and (k + q) 2 = z, we have that k ■ q = ^Jxyco^x-, and so 
(k ■ q) 2 j q 2 = y cos 2 x, and z = x + y + 2^/xycos \- Due to the dependence of the unknown 
function on A(z), to solve (13.31) one should carry out (numerically) the two integrals defined 
in (13.41) . Alternatively, one resorts to standard approximations (see ( 1A-4I) ) for the integral 
over dx, thus reducing the numerics to a single one-dimensional integral. In our analysis we 
will use a modified version of the usual angular approximation (see Appendix). 

The first important issue is the value that (13. 3p furnishes for A _1 (0); in particular, in 
order to obtain IR-finite solutions, A _1 (0) should be non- vanishing. It turns out that A _1 (0) 
may be determined from (13.31) exactly, before doing any approximation. Specifically, 



d q k- 



-A(k)A(k + q) 



2tt / dx sin x cos \ \ dyy A (y) 



l - J d A kk 2 A 2 (k). 



Thus, one obtains from (13.31) 

36<? 2 



A" l (0) 



5tt 2 



2(1 + ci) J d 4 kA(k) - (1 + 2ci) J d A kk 2 A 2 (k) 



(3.5) 



(3.6) 
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Since the r.h.s. of (13. 6p is divergent, in the next section this expression will be appropriately 
regularized following the rules of dimensional regularization; there it will become clear 
why the UV running of the effective mass is of central importance. 

Having determined the exact expression A _1 (0), one might be tempted to carry out the 
angular integration by employing the usual approximation of (IA-40 ; however, particular care 
is needed, since the straightforward application of (1A-41) would misrepresent quantitatively 
some essential features of the original Eq. fl3.3p . and most importantly the terms determin- 
ing the running of the mass. To remedy this, in the Appendix we implement judicious 
modifications to some of the relevant results of the angular approximation. 

Let us denote the r.h.s. of (13. 3p by I(q 2 )- Then we first write (13.31) schematically in the 
form 

A _1 (x) = [I{x)-I{0)]+A- 1 {0) 

« [I MA (x)-I MA (0)] + A-\0), (3.7) 

where the exact A _1 (0) is given by (13.61) . and the subscript "MA" denotes that the (modified) 
angular approximation has been employed. Then, using for all terms appearing in ( I3.3P the 
expressions given in ( 1A-3h and ( 1A-21I) . after some algebra we arrive at the SDE 



A-^ar) = Kx + bg 2 ^A i {x) + A-\0) 



(3.8) 



i=i 



with 



A x {x 
A 2 (x 
A 3 (x 
A 4 (x 
A 5 (x 
A e (x 
A 7 (x 
A 8 (x 



6c 2 



POO 

x / dyyA 2 (y) 

J X 



^x J dyA(y), 

-~(l + c 1 jA(x) J*dyy 2 A(y), 

X 

dyyA(y) , 

X 

dyy 3 A(y), 



Qc 2 
2 A(x 



5 x 



i [ X dyy 3 A 2 (y) 



(3.9) 
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and K is the wave-function renormalization constant, whose value is fixed by the renormal- 
ization condition A _1 (/i 2 ) = /z 2 , with /j, 2 3> A 2 , where A is the mass-scale of QCD. 

As has been discussed extensively in the literature, due to the Abelian Wis satisfied 
by the PT effective Green's functions, A _1 (g 2 ) absorbs all the RG-logs 0, H, Q, 0, 



exactly as happens in QED with the photon self-energy. Equivalently, since Z g and Z^, 
the renormalization constants of the gauge-coupling and the gluon self-energy, respectively, 



-1/2 



17|, the product d(q 2 ) = g 2 (fj, 2 )A(q 2 , /i 2 ) forms a 



satisfy the QED relation Z g = Z~ 
RG-invariant (//-independent) quantity; for large momenta q 2 , 

d(q 2 ) = 92pert [ q2) , (3.10) 
q 2 

where g^ ert {q 2 ) is the perturbative limit of the RG-invariant effective charge of QCD, i.e. 

2 / 2n _ 9 2 {^) _ 1 mn 

9 P erAQ ) i + hg 2 {^ 2 )\n{q 2 /^ 2 ) Mn(g 2 /A 2 )' l ' ' 

Notice however that Eq. (13.81) does not encode the correct RG behavior, in the sense that 
it cannot be cast in a form containing only RG invariant quantities, such as d(q 2 ) and 
Qpertil 2 ) ■ This can be ultimately traced back to the fact that the various gauge technique 
inspired Ansatze for the full vertices (in our case the three-gluon vertex) tend to mishandle 
their transverse (identically conserved parts); this, in turn, forces one to renormalize sub- 
tractively instead of multiplicatively. We emphasize that there exist systematic methods 
for the construction of transverse parts with the correct UV properties [3^], but have not 
been extended to the case of the three-gluon vertex. Therefore, to restore the correct RGI 
properties, we follow the heuristic procedure first proposed in [2J, [25] , used recently also in 



Specifically, every A(t) (t = y,x) appearing inside the integrals on the r.h.s. of (13. 8p is 
to be multiplied (by hand) by a factor 1 + bg 2 ln(t / fi 2 ) = g 2 /gl ert (t), i.e. we carry out the 
replacement A(t) — > [g 2 / gp ert (t)~\ A(t) (with b — ► b). Then, one may rewrite (I3.8P in terms 
of two RG invariant quantities, d(q 2 ), and 

C(q 2 )= 9p - 2 t (q 2 )=b\n(q 2 /A 2 ) , (3.12) 

as follows: 

8 

d-\x) = K'x+bJ2M x ) +d~ 1 (0), (3.13) 

i=l 
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with 



AAx) 



l + ^)x I dyyC 2 {y)d\y) 



6co 



A 2 (x) = —x I dyC(y)d(y) 



A 3 (x) 
AAx) 



6c \ f x 
1 + xC{x)d{x) J dyyC(y)d{y), 



6 



A 5 (x) = --^l + Cl jC(x)d(x) J^dyy 2 C(y)d(y) 



A 6 {x) 



Qc 2 

~ Jo 



dyyC(y)d(y) 



A 7 (x) = 2 -C{x) d -^- rdyy 3 C(y)d(y), 

X Jq 

M*) = ir [ X dyy 3 £\y)d 2 (y), 



5.r 



where K' = K/g 2 , given in closed form by 



K' = l 



8=1 



and 



with 



d-\o) 



3b 

57T 2 



To 



2(1 + ci)T - (1 + 2c x )T 1 

d 4 k C(k 2 ) d(k 2 ) , 
d A kk 2 C 2 (k 2 )d 2 (k 2 ). 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



IV. THE UV BEHAVIOR OF THE EFFECTIVE GLUON MASS 

In this section we will discuss in detail the UV behavior of the effective gluon mass 
obtained from the SDE of f)3.13p . As we will see, depending on the values of the parameters 
c\ and C2, originally appearing in the vertex Ansatz of (13. 2p . one obtains logarithmic or 
power-law running as two distinct dynamical possibilities. 

We begin by briefly reviewing the importance of the running of the gluon mass for ob- 
taining a finite value for d -1 (0); the finiteness of <i -1 (0) is intimately linked to the renor- 
malizability of the theory; if <i _1 (0) turned out to be divergent, there would be no consistent 
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way to eliminate this divergence by absorbing it into the renormalization of the parameters 
appearing in the fundamental QCD Lagrangian. 

As has been explained in 8j, under special assumptions on the form of d(k 2 ), the r.h.s. of 
(I3.16P can be made finite by simply employing standard dimensional regularization results. 
Specifically, let us write d(q 2 ) as 



d(q 2 ) = <?V)A(^) 



where 



and 



g 2 ^ 



A(g 2 



b In 



1 



q 2 + m 2 (q 2 ) 

q 2 + f(q 2 ,m 2 (q 2 )) 
A 2 



-i 



(4.1) 



(4.2) 



(4.3) 



For large values of q 2 , A(q 2 ) — > 1/q 2 and g 2 (q 2 ) — > gp ert (<l 2 ), an d d(q 2 ) reduces to its 
perturbative expression of ( I3.10p . Then, it is relatively straightforward to show, using the 
elementary result 



dz 



1 



z (lnz) 1+ T 
that the difference T — T 1 , given by 



7 (In z)' 



(4.4) 



d 4 kC 2 (k 2 )m 2 (k 2 )d 2 (k 2 ) 



+ b I d A k £(k 2 ) g 2 {k 2 ) d{k 2 ) In ( 1 + 



f(k 2 ,m 2 (k 2 )) 
k 2 



(4.5) 



is finite, provided that m 2 (k 2 ) drops asymptotically at least as fast as ln~ a (k 2 ), with a > 1, 
and f(k 2 ,m 2 (k 2 )) as ln~ c (k 2 ), with c> 0. 

Then, T may be regularized simply by subtracting from it J d A k/k 2 = 0, i.e. 

1 



1 



cffc C(k 2 )d(k 2 ) - 



h m 2 (k 2 )A(k 2 ) 



k 2 



d^k 



k 2 



b / d*k d{k 2 ) In 1 



f(k 2 ,m 2 (k 2 )) 
k 2 



(4.6) 



Thus, the regularized expression for d l (0) is given by 

36 



<(0) 



5tt 2 



Tl eg + (l + 2c l )(T Q -T 1 ) 



(4.7) 



When solving (13. 13[) one must impose (14.71) as an additional constraint; we will refer to (14. 7p 
as the "seagull-condition" . The way this is done is by first choosing an arbitrary finite value 
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for the c/ _1 (0) appearing in (13 . 13[) and then solving it numerically. The solution obtained for 
d(q 2 ) must be first decomposed following (14. II) - (14.3 1) , and then be substituted into the r.h.s. 
of (14. 7p ; self-consistency requires that the resulting expression for <i~g(0) must coincide with 
the initial value d~ 1 (0). Evidently, the way that m 2 (q 2 ) runs is essential for this procedure, 
and can affect considerably the quantitative predictions. 

In order to obtain from (13.131) the equation that determines the behavior of m 2 (x) at large 
x, first set in the r.h.s. of (I3.14p xC(x)d(x) — > 1, C(x)d(x) — > 1/x, and C{y)d{y) = A(y). 
Next, use the identity yA(y) = 1 — m 2 (y)A(y) in all Ai(x), keeping only terms linear in 
m 2 (terms quadratic in m 2 are sub leading and may be safely neglected). Then separate 
all contributions that go like x from those that go like m 2 on both sides, and match them 
up. This gives rise to two independent equations, one for the "kinetic" term, which simply 
reproduces the asymptotic behavior x In x on both sides, and one for the terms with m 2 (x), 
given by 

f x 0,2 f x 
m 2 (x)\nx = 6 _1 rf _1 (0) + a\ \ dy m 2 (y) A(y) H / dyym 2 (y)A(y) 

Jo x Jo 

+ ^[ dyy 2 m 2 (y)A(y) + a i x [ dym 2 {y)A 2 {y) , (4.8) 



o 



with 



(3 . 4 6ci 2 6c 2 > 

a x = -(1 +c 2 - d) , a 2 = - + — , a 3 = --, a 4 = 1 + — . (4.9) 



Then, rewrite the first integral on the r.h.s of (14.101) as JJf = J °° — to obtain 

f°° a<i f x 
m 2 (x)lnx = C — ai / dym 2 (y)A(y) H / dyym 2 (y)A(y) 

Jx x J0 

+^ [ dyy 2 m 2 (y)A(y) + a i x [ dym 2 {y)A 2 {y) , (4.10) 

x JO Jx 



with 



C = r i rf- i (0) + a 1 / dym 2 (y)A(y). (4.11) 



o 



As we will see in a moment, the two possible asymptotic solutions of physical interest for 
m 2 (x) are given by 

m\{x) = A?(lnx)- (1+7l) , (4.12) 
ml{x) = ^(lna;) 72 - 1 , (4.13) 

where Ai and A 2 are two mass-scales, and ji > 0, i = 1,2. The possibility of power-law 
running for the effective gluon mass, as expressed by (14.131) . was first conjectured in j^, 
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motivated by similar results in the study of chiral symmetry breaking. Indeed, notice the 
similarity between the solutions given in (I4.12p and (14.131) for the effective gluon mass and 
those appearing in the more familiar context of the SDE (gap equation) for the quark self- 
energy; there, one finds the following two asymptotic solutions for the dynamically generated 
quark mass M(x): 

Mi(x) = /xi(lnx)" 7 ' , 

M 2 (x) = ^(lnx) 7 '- 1 , (4.14) 

x 

where 7/ = 3C//167r 2 6, with Cf the Casimir eigenvalue of the fundamental representation 



[C f = (N 2 - 1)/2N for SU(N)] [38]. 



Now, the important point to appreciate is that, in order for ( 14.101) to have solutions 
vanishing in the UV, it is necessary that the constant term on the r.h.s. vanishes, i.e. C = 0. 
If, for some reason, this condition cannot be implemented, the solution obtained will reach a 
constant value in the deep UV, thus invalidating the basic characteristic of the dynamically 
generated mass. Given that both c? _1 (0) and the integral appearing in (14. lip are manifestly 
positive quantities, one obvious necessary condition for obtaining C = is that a x < 0. This 
requirement, in turn, restricts the possible values of the parameters c\ and C2, through the 
equation defining a\ (first in ( I4.9P ). Assuming the correct sign for a\, the way to actually 
enforce C = will be completely dynamical: one must look for masses with the appropriate 
momentum dependence such that the r.h.s. of (14.111) can be made equal to zero. As we will 
see in the next section, the condition (14. lip [or, its improved version, (I4.28p ] constrains the 
behavior of the dynamical mass in the IR and intermediate momentum regimes. 

Next, we set C = in ( I4.10p and verify that indeed mi(i) or m 2 (x) [Eqs.f l4.12~]) and ( 14.131) ] 
satisfy it. The upshot of this analysis will be that when mi (2) is substituted into the r.h.s 
of (14.101) the first integral provides the solution, while all others are subleading, whereas 
for 777,2(0;) the leading contribution comes from the second integral, and the other three are 
subleading. (The third and fourth integrals are thus subleading for both types of solutions). 

In our demonstration we will employ the asymptotic property of the incomplete T func- 



tion. The latter is defined as 39 ] 



00 

t ,a-l 



T(a,u) = / dte-'e-' , (4.15) 

J u 

(with no restriction on the sign of a), and its asymptotic representation for large values of 
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I it I is given by 

T{a,u) = u a - 1 e- u + 0{\u\- 1 ). (4.16) 

To see in detail what happens in the case of the logarithmic running, substitute (14.1 2p 
into both sides of (I4.10p and use that asymptotically A(y) — > y^ 1 . Then 

dy = 7 r x m\ (x) In x , (4.17) 



X 

X 



dym\(y) = m\{x) + O (l/lnz) , (4.18) 

- 2 [ X dyyml(y) = ™M + (D (l/]nx) , (4.19) 
x Jo * 

x / dy-^- = m\{x) + 0(ll\nx) . (4.20) 

Jx y 2 

In evaluating the integral of ( 14. 17ft we have used (14.41) ; for the remaining three we have set 
y = e~ l in ( 14. 18ft . y = e _t//2 in ( 14. 19(1 . and y = e* in ( 14.20(1 . to cast them into the form 
of the incomplete F function given in (14.151) . and have subsequently used the asymptotic 
expression of (14.161) . 

Thus, for asymptotic values of x the dominant contribution comes from (14.171) Substitut- 
ing into (14.101) . we see that both sides can be made equal provided that 

7i = -a t . (4.21) 

Since 71 must be positive, it follows that a\ < 0, consistent with the requirement imposed 
by Eq. ( 14.111) . as discussed after Eq. ( 14.14(1 . 

Let us next turn to the case of power running, substituting (14. 13(1 into both sides of 
(14. 10(1 . Now the leading contribution comes from the integral proportional to 02 in (14.101) . 
Specifically, 

f°° dyH^M. = m 2 2 (x) + 0(l/\nx) , (4.22) 

j x y 

- / dymliy) = l2 lm l( x ) lna; + - > ( 4 - 23 ) 
x l n ~ x 



X 2 



dyym\{y) = ml{x) + O (l/lnx) , (4.24) 



00 



x I dy^fi- = *l + 0(l/lnx) . (4.25) 

Jx y 2 

In evaluating (14.231) we have used (I4.4p . The constant d comes from the lower limit of 
the integral; it is finite, because in that limit one must use inside the integral the full 

17 



yA(y), which is infrared safe due to the presence of the mass. The term proportional to 
d is suppressed by a factor In 72 x (assuming 72 > 0) compared to the first term, and can 
therefore be neglected. As in the case of the logarithmic running, for the other three integrals 
we have used the appropriate change of variables to cast them into the incomplete V function, 
resorting again to its asymptotic expression. Thus, we conclude that m 2 (x) satisfies (I4.10p 
provided that 

72 = a 2 . (4.26) 

This condition, in turn, constrains the possible values of c\ appearing in the definition of a 2 
(second of (Oil ). 

We emphasize again that, for either of the two physically relevant possibilities given by 
(14. 1211 and (14.131) to be realized, the constant term on the l.h.s of (14. Ill) must be forced to 
vanish, by imposing the mass-condition 

POO 

d-\0) = ll b dym 2 (y)A(y). (4.27) 







In addition, it is important to mention that whereas the individual ( 14. 121) and ( 14.131) are 
separately solutions of (14.101) . a linear combination of the form m 2 (x) = Cim\{x) + C^mf (x) 
cannot be regarded as a solution. The reason is that for m\{x) to be a solution one neglects 
all terms in (I4.18I) - (I4.20I) . which are, however, clearly larger than m\{x). This does not 
necessarily mean that the two runnings cannot coexist, it simply says that the possible 
coexistence cannot be self-consistently inferred from (14.101) . For this reason, in the analysis 
of the next section the two possibilities will be treated separately. It should also be clear 
that the terms quadratic in m 2 (x) that have been dropped when deriving (14.81) are indeed 
subleading for both types of asymptotic behavior, (14.121) and ( 14. 13D . 

Since in this section we have been mainly interested in the UV running of m 2 (x), in the 
analysis presented above we have used for the g 2 (q 2 ) of (14. 3[) its perturbative expression, 
given in ( 13. lip , i.e. we have replaced g 2 (q 2 ) — > glertiQ 2 )- A complete treatment, where the 
full expression for g 2 (q 2 ) is kept, does not affect in the least the conclusions regarding the 
UV behavior of m 2 (x), but modifies slightly the condition ( I4.2T[) in the IR. Specifically, 



00 



d-\0) = (l + ll )~b / dyym 2 (y)£ 2 (y)d 2 (y)-~b / dym 2 (y) C(y)d(y) . (4.28) 
Jo Jo 

Evidently, in the limit C{y) d{y) — > A(y), (I4.28P reduces to (14.271) . as it should. In the numer- 
ical analysis that follows we will always use (I4.28[) . referring to it as the "mass- condit ion" . 
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V. NUMERICAL ANALYSIS 



In this section we will solve numerically the integral equation given in Eq. (13. 131) supple- 
mented by the renormalization condition (13.151) . and subject to the two constrains imposed 
by the mass- and seagull-conditions, Eqs. (14.281) and (14.71) . respectively. As mentioned already 
in the previous sections, and as we will see in detail in what follows, the first condition re- 
stricts the momentum-dependence of the mass in the intermediate and deep infrared regimes, 
while the latter furnishes essentially the value of <i _1 (0). 

When dealing with this problem we have at our disposal three undetermined parameters: 
Ci and C2 appearing in the Ansatz for the three-gluon vertex [Eqs. (13.11) and (13. 2p ], and the 
value of <i _1 (0). It turns out that the simultaneous solution of the integral equation and its 
constrains restricts considerably the acceptable combinations of these parameters. 

The strategy we will employ in our numerical analysis consists of the following main 
steps: 

(i) We choose an arbitrary initial value for d -1 (0), to be denoted by 0^(0), together 
with a set of values for ci, c 2 , and we substitute them into the integral equation (13.131) . 
generating a solution for d(q 2 ). 

(ii) The solution for d(q 2 ) obtained in (i) must be then decomposed as the product of 
A(g 2 ) and g 2 {q 2 ), according to Eqs. (l4.ip . (14.21) . and (14. 3p . To do this, first a simple Ansatz 
for m 2 (q 2 ) is written down, which in the UV displays one of the two physically relevant 
asymptotic behaviors (logarithmic or power-law), while in the IR reaches a finite value. The 
generic form of these two types of Ansatze is given in (15.11) and (15. 3p . for logarithmic and 
power-law running, respectively. The anomalous dimensions 71 and 72 appearing there are 
linear combinations of c\ and C2 [given by (I4.9p . (l4.2ip . and (14.261) ]. and control the behavior 
of m 2 (q 2 ) in the deep UV, whereas the parameters p and mo are free for the moment, and 
affect the intermediate and IR regions. 

(iii) The integrals on the r.h.s. of the mass-condition (14.281) are evaluated numerically, 
using as input the Ansatz for m 2 (q 2 ) chosen in (ii), together with the numerical solution for 
d(q 2 ); this furnishes a value for the g? -1 (0) on the l.h.s. By varying p and m we try to make 
that g? -1 (0) match £^(0) of step (i). We may or may not be able to do this, depending 
on the values of C\ and C2, and the type of Ansatz (logarithmic or power-law) chosen for 
the mass. In general, in the cases where (14. 28[) can be satisfied, one finds that this may be 
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accomplished not just for one but for various sets of values for p and m®. This, in turn, 
gives rise to a family of possible masses, which have a common behavior in the deep UV, 
but differ in the intermediate and IR regions. 

(iv) Once the family of allowed masses has been determined from the mass- condition in 
(iii), one extracts from (14.11) the corresponding families of non-perturbative effective charges 
g 2 {q 2 )', one simply multiplies the numerical points of d(q 2 ) by [q 2 + m 2 (q 2 )}. On physical 
grounds we require that the resulting effective charges should be monotonically decreasing 
functions of q 2 ; this requirement eliminates any member of the mass family that gives rise 
to effective charges with "bumps". All g 2 (q 2 ) so obtained display in the UV the logarith- 
mic running expected from the one-loop RG (i.e. asymptotic freedom corresponding to 
(3 = —bg 3 ), and reach a finite value in the deep IR. For the implementation of the seagull- 
condition, (14. 7p . one must supply the functions / (q 2 , m 2 (q 2 )) , appearing in (14.31) . Therefore, 
we must extract for each g 2 (q 2 ) the corresponding f{q 2 ,fn 2 (q 2 )); the way this is done is by 
fitting the numerical points determining g 2 (q 2 ) by Eq. (14.31) . assuming for f(q 2 ,m 2 (q 2 )) the 
expression of (15.21) . 

(v) We next substitute g 2 (q 2 ), m 2 (q 2 ) and d(q 2 ) into the integrals on the r.h.s. of 
the seagull-condition (14. 7p . whose value is computed numerically; this furnishes the value 
for <i _1 (0) appearing on the r.h.s. If this value for <i _1 (0) does not coincide with d~^(0) [we 
require an accuracy of about 1 part in 10 3 ] a new set of values for c\ and C2 is chosen (keeping 
d~^(0) fixed), and the procedure is repeated from the beginning, until coincidence has been 
reached. At that point we consider to have found a solution, namely the d(q 2 ) obtained for 
din{0) and the values for c\ and C2 used the last (and only "successful") iteration. 

(vi) A different value for d~^(0) is chosen, and the procedure is repeated starting from 
step (i). 

To solve the integral equation we employ a simple iterative procedure, where an initial 
guess is made for the solution d(q 2 ) on a discretized momentum grid in the domain of [0, Auv] • 
More specifically, the grid is split in two regions [0, /i 2 ] and (/z 2 , Auv] whose purpose is allow 
for the implementation of the renormalization condition, given by Eq. ( 13. 1 5f) . Typically, 
we choose Auv — 10 6 GeV 2 , /j 2 = M| = (91. 18) 2 GeV 2 and we used as input a value of 
A = 300 MeV for the QCD mass-scale. 

Our numerical analysis reveals a clear separation between the two types of asymptotic be- 
havior for m 2 (q 2 ) depending on the values chosen for c\ and C2. Specifically, for c\ G [0.15, 0.4] 
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and C2 G [—1.07, —0.92] the asymptotic behavior of m 2 (q 2 ) is given by Eq. (l4.12l) . whereas 
for ci G [0.7,1.3] and c 2 G [—1.35,-0.68] the m 2 (q 2 ) displays the power-law running of 

Eq.flHSD. 

A. m 2 (q 2 ) with logarithmic running 

When solving Eq. (13.131) choosing values for c\ and c 2 from the intervals C\ G [0.15,0.4] 
and C2 G [—1.07,-0.92], the constraints (14.281) and (14.71) can be simultaneously satisfied, 
and the d(q 2 ) obtained may indeed be decomposed as in Eq. (l4.1l) . with a functional Ansatz 
for the running mass of the form 



2i 2\ 2 

m (q ) = m 



ln , q 2 + pm 2 \ a ( pm 2 



A 2 I V A 2 



-I-71 



(5.1) 



where 71 = —d\ [see (14.211) and (14.91) ]. Evidently, for large q 2 the above expression goes 

r / 2 \ "1 1+71 

over to the logarithmic behavior described by Eq. (14.121) , with A 2 = m\ ln I 



This 

simple Ansatz connects continuously the UV and IR regions; at q 1 = reaches the finite 



x A2 

2 



m, 



value m 2 (0) 

The parameters mo and p appearing in Eq. (15.11) control the way the mass runs in the 
intermediate and IR regions; their values are restricted by the mass-condition, Eq. (14.281) . 
To impose the mass-condition, we first choose a random value for m , and then we search 
for values of p that satisfy Eq. (14.281) . Even though this procedure does not single out 
a unique pair of values for mo and p, it restricts considerably their allowed range. In 
fact, the acceptable range for (mo, p) gets further restricted by imposing the additional 
requirements that all the m 2 (q 2 ) and the effective charges generated subsequently from 
them (by multiplying d(q 2 ) by [q 2 + m 2 (q 2 )]) should be monotonically decreasing functions 
of q 2 . The combination of all these constraints leads eventually to rather stable results: if a 
pair (m ,p) furnishes a consistent solution for a given (01,02) , then any other pair (m' Q , p') 
is also a solution, provided that C2 is only slightly adjusted (less than 5%). 

The running couplings g 2 (q 2 ) obtained using Eq. (l5.ip in Eq. (14.11) can be fitted very ac- 
curately by means of Eq. (l4.3p . with the function / '(q 2 , m 2 (q 2 )) given by 

f(q 2 , m 2 {q 2 )) = Pl m 2 (q 2 ) + p 2 + P 3 ]2 f ( f ( \^ , (5-2) 

q z + m z {q z ) [q z + m z {q z )\ A 

In Fig.(T5]) we present a typical solution for d(q 2 ), m 2 (q 2 ), and the effective charge 
a (l 2 ) — 5 ,2 ( ( ? 2 )/4vr, respectively. 
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q [GeV] 





q 2 [GeV] 

FIG. 2: Solutions obtained for the choice d -1 (0) = 0.02 GeV 2 , c x = 0.15 and c 2 = -0.9635. Upper 
panel: the numerical solution for d(q 2 ). Lower panels: the logarithmic dynamical mass, m 2 (q 2 ), 
for uiq = 0.3 GeV 2 and p = 1.007 in the Eq. (|5.ip . On the right panel we show the running charge, 
a(q 2 ) = g 2 (q 2 )/4ir, which can be fitted using Eqs.flOJ) and §52§ with p x = 6.378, p 2 = -8.984 and 
p 3 = 3.466. 



B. m 2 {q 2 ) with power-law running 

As we increase c\, Eq. (14.281) can not be satisfied if we insist on imposing the loga- 
rithmic running for m 2 (q 2 ). In particular, for c\ G [0.7,1.3], c 2 G [—1.35,-0.68] and 
d-\0) G [0.01 GeV 2 , 0.04 GeV 2 ] we have verified that the mass-condition can be satisfied 
only if instead of Eq. fl5.ip we use the following functional form for the running mass 

72-1 



m 2 (q 2 



rrr n 



q 2 + rriQ 



In 



q 2 + p ml 
A 2 



In 



A 2 



(5.3) 



with 72 = ci2 [see (14.261) and (14.91) ]. Evidently, this Ansatz corresponds to power- law running 

for the mass; for large q 2 , the m 2 (q 2 ) goes over to the solution denoted by m\(q 2 ) in (14.131) . 

1 1-72 

. Clearly, (15. 3p is the simplest extension of the asymptotic 



with A2 



In 



expression (I4.13P to the entire range of momenta, from the UV all the way down to q 2 = 0, 
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where it assumes the finite value m 2 (0) = m 2 ,. 

As in the previous case, only those sets of mo and p that satisfy the mass-condition (I4.28P 
are allowed. In fact, in the case of the power-law running this condition turns out to be 
significantly more restrictive than in the logarithmic case. The reason is that, due to the 
faster decrease of the mass in the UV, the leading contribution to the mass-condition comes 
now from the intermediate and IR regions; therefore, the result is much more sensitivity to 
small variations of mo and p. 

A typical solution for a choice of ci, c 2 , and c? _1 (0) within the aforementioned ranges is 
shown in Fig. ([3]). The d(q 2 ) is decomposed according to ( 14. II) into an m 2 (q 2 ) of the general 
form given in Eq. fj5.3j) and an effective charge a(<? 2 ); the latter is fitted using again (14. 3 j) 
and Q . 




Gluon Propagator 



10"* 10"' 10 1 10 J 

q 2 [GeV 2 ] 




10" 2 10"' 10° 10' 10 2 10"' 10"' 10' 10 3 



q 2 [GeV 2 ] q 2 [GeV 2 ] 

FIG. 3: Solutions obtained for the choice c2 _1 (°) = °- 04 GeV 2 , c x = 1.1 and c 2 = -1.121. Upper 
panel: the numerical solution for d(q 2 ). Lower panels: the power-law dynamical mass, m 2 (q 2 ) for 
to 2 , = 0.5 GeV 2 and p = 1.046 in the Eq. (|5.3p . On the right panel we show the running charge, 
a{q 2 ) = g 2 (q 2 )/Air, which can be fitted by Eqs.jOJ and with Pi = 1-205, p 2 = -0.690 and 
p 3 = 0.121. 



23 



In the upper panels of Fig.Q we plot a series of effective charges obtained by fixing 
different set of values for rf _1 (0), ci, and Ci- All these coupling were subjected to the 
constraints imposed by Eqs. (14.281) and (14. 7ft ; the corresponding values for C\ and c 2 are 
given in the legend. The respective logarithmic and power-law masses are shown in the 
lower panels. Observe that a(0) shows a strong dependence on the values chosen for c\ and 
C2~, this is so, even if the same value for <i _1 (0) is chosen. The difference between effective 
charges obtained with the same (i _1 (0) is due to the fact that one is forced to use different 
values of in order to satisfy the mass-condition. Thus, by changing the value of Mq, one 
obtains different values of a(0) for the same (i _1 (0). 




Running Charge - Logarithmic mass 
-♦-d 1 (0)=0.03, c,=0.40 and c 2 = -0.988 
-'-d" 1 (0)=0.04, c^O.30 and c 2 = -0.995 
-— d 1 (0)=0.02, c,=0.20 and c 2 = -0.964 
— ■- d" 1 (0)=0.05, c =0.25 and c = -1.014 
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Logarithmic running mass 

-.— d' 1 (0)=0.03, ^=0.40 and c 2 = -0: 
-•- d 1 (0)=0.04, c,=0.30 and c = -0.995 
-•-d' , (0)=0.02, c^O.20 and c 2 = -0.964 
— d"'(0)=0.05, c, =0.25 and a = -1.014 





Running Charge - Power law mass 

d"'(0) = 0.04, d = 0.7 and c2 = -0.838 
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FIG. 4: Upper panels: The running charges, a(q 2 ) = g 2 (q 2 ) / Air , corresponding to the different 
choices of (i _1 (0), c\ and c 2 for the logarithmic running mass (left panel) and for the power-law 
running mass (right panel). The corresponding masses are plotted in the lower panels. 
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In Fig.© we show a comparison between logarithmic and a power-law running, by choos- 
ing ci and C2 from the corresponding intervals, for fixed values of <i _1 (0) and tUq. The corre- 
sponding d(q 2 ) are shown in the upper panel, whereas the running masses and couplings are 
plotted in the lower left and right panels, respectively. The faster decrease of the power-law 
running mass in the UV is clearly visible. It is evident that, as already mentioned, in the 
power-law case the leading contribution to the mass-condition comes mainly from the IR and 
intermediate regions, while in the logarithmic case the UV region provides a considerable 
support. Since the UV behavior of the two d(q 2 ) and a(q 2 ) is essentially fixed by asymptotic 
freedom, whereas their IR regimes are determined, to a large extent, by the value of tuq, the 
difference between the two cases is perceptible only in the intermediate momentum region. 




10' 3 10'' 10' 10 3 10' 3 10"' 10' 10 s 



q 2 [GeV 2 ] q 2 [GeV 2 ] 

FIG. 5: Comparison between logarithmic and a power-law running cases. Upper Panel: The nu- 
merical solutions for d(q 2 ), when d~ 1 (0) = 0.02, GeV 2 , c\ = 0.15 and C2 = —0.9635 (Logarithmic 
case) and for c\ = 1.00 and C2 = —1.12 (Power-law). Lower Panels: The corresponding logarithmic 
and a power-law running masses and the running charges. 
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VI. DISCUSSION AND CONCLUSIONS 



In this article we have studied both analytically and numerically the non-linear SDE for 
the gluon self-energy in the PT-BFM formalism, focusing on the dynamical generation of 
a gauge-invariant infrared cutoff. In particular, we have established the existence of IR- 
finite solutions, i.e. solutions that are finite in the entire range of momenta, displaying 
asymptotic freedom in the UV and reaching a finite positive value at q 2 = 0. This non- 
perturbative behavior may be described in terms of an effective gluon mass, whose presence 
tames the perturbative Landau singularity and prevents A(g 2 ) from diverging in the IR. 
Just as happens with the constituent quark masses, this dynamical gluon mass depends 
non-trivially on the momentum. Our study of the non-linear SDE has revealed that the 
dynamical mass m 2 (q 2 ) may have two different types of functional dependence on q 2 in the 
deep UV: (i) m 2 (q 2 ) drops as an inverse power of a logarithmic; this behavior has also been 
found in the studies of linearized SDE, and (ii) m 2 (q 2 ) with power- law running, i.e. the 
mass drops off as 1/q 2 . This type of solution is found for the first time in the context of a 
SDE. At the level of the SDE we study either type of asymptotic behavior (logarithmic and 
power-law) may be obtained, depending on the details of the three-gluon vertex. The latter 
depends on two parameters, c\ and C2, which control the relative contribution of its various 
tensorial structures. Our numerical analysis reveals that the sets of values for c\ and C2 that 
give rise to logarithmic running belong to an interval that is disjoint and well-separated from 
that producing power-law running. 

The possibility of gluon masses falling like the inverse square of the momentum has been 
anticipated Q by analogy with the constituent quark masses generated from the standard 



gap equation for the quark self-energy [second equation in f)4.14p ]. In addition, general OPE 
considerations support the existence of a m 2 (q 2 ) displaying power-law running. Assuming 
that the OPE holds for a quantity like m 2 (q 2 ), and given that, in the absence of quarks, 
(G 2 ) is the lowest order (dimension four) local gauge-invariant condensate, then one would 
expect that asymptotically, and up to logarithms, m 2 (q 2 ) ~ {G 2 )/q 2 , exactly as was found 



To be sure, the various connections between the effective gluon mass and the OPE, the 
gluon condensates, and the QCD sum rules, deserve a detailed, in-depth study. One issue 
is the type of modifications induced to the OPE predictions for observables (i.e. correlators 



in [26]. 
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of gauge-invariant currents) if one were to use in their calculation gluon propagators with 
a dynamical (or even hard) mass, as was first done in 40j. In addition, despite important 



contributions in this directio n |26l . |41|. |42||. a definite, first-principle relation between m 2 and 
(G 2 ) (or other condensates 45j) still eludes us. In the context of the SDE this is mainly 
because the CJT formalism [43j has not been yet fully adapted to treat gluon mass generation 
in a consistent way. Qualitatively speaking, the CJT effective potential V is given by 

V = ~Tr In (AAo 1 ) + l - (TrAA 1 - l) + V 2PI (6.1) 

where the trace is taken in the functional sense, and V<jpi denotes the contributions from the 
(appropriately dressed) two-particle irreducible graphs. In the original formulation V is a 
functional of the conventional gluon propagators, and higher point Green functions; its ex- 
tremization with respect to any of yields the corresponding SDE's. To make reliable contact 
with the results of the BFM-PT, one should modify Eq.f l6.ll) appropriately, expressing it in 
terms of the gauge-invariant PT gluon propagator. Thus, the gluon mass will enter into V 
through the massive gluon propagators; then, the minimization of V will yield a theoretical 
expression for the energy density of the QCD vacuum, which must be set equal to the exper- 
imental value obtained using QCD sum rules. To date, the aforementioned modifications to 
V have been carried out at the two- loop level only ; clearly, their all-order generalization 
would be of great interest. 

It is clear from the analysis presented that the role of the three-gluon vertex vertex is 
absolutely central. Specifically, the tensorial structure of the vertex, the presence or absence 
of kinematic poles, and the relative strength between the various components are determining 
factors for the existence of IR-finite solutions, and the type of running of the effective gluon 
mass. As we have mentioned in Sec. the Ansatz of (13.11) employed for the vertex attempts 
to capture some of the main features, but should be eventually obtained from an independent 
study of the dynamical SDE that it satisfies. The rich tensorial structure of a vertex with 
three Lorentz indices turns such a study into a rather complicated task, given that, in 
general, one has to deal with fourteen form-factors. However, as a first approximation, 
one can focus on those form-factors that enter into the expression determining A _1 (0), i.e. 
the generalized version of (13.61) . As an alternative, one could try to improve the Ansatze 
employed for the vertex, in the spirit of [371 ]. in an attempt to correctly incorporate the 
required asymptotic behavior into the SDE, i.e. without having to resort to the heuristic 
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procedure followed here. 

Another source of relative uncertainty when evaluating the SDE in question is the use 
of the angular approximation. As we explained in detail, the standard version of this ap- 
proximation gives rise to an approximate SDE that does not capture faithfully some of 
the essential features of the original SDE. To ameliorate this drawback we have introduced 
modifications to the angular approximation, presented in the Appendix. In our opinion the 
qualitative conclusions of this article are robust and do not depend on its use; we expect 
them to persist a more complete study, where the angular integration is carried out numeri- 
cally, without resorting to any approximation. Should quantitative discrepancies arise, they 
will mainly manifest themselves in changes of the values of c\ and C2, which, at this level of 
approximation, are free to vary. To be sure, once dynamical information for the three-gluon 
vertex has been furnished one would be more restricted, and the two-dimensional integration 
should be carried out in its entirety. 

Last but not least, let us turn to the possible role of the ghost sector. As we have amply 
emphasized, in the PT-BFM scheme the omission of the ghost loops does not interfere with 
the gauge- invariance of the final answer, and in particular with the transversality of the gluon 
self-energy. In addition, neglecting ghost contributions only affects the RG-logarithms by 
about only 10%, the difference between b = 11Ca/487t 2 and b = IOCa/^tt 2 . The possible 
impact of the ghosts in the IR is, however, an entirely different matter; it will greatly depend, 
among other things, on the structure and solutions of the corresponding SDE for the ghost 
propagator within the PT-BFM formalism. Actually, given that in this latter framework one 
is working in the Feynman gauge (of the BFM) there is no a-priori reason that would exclude 
the possibility of obtaining IR-finite solutions for the ghost propagator. If this turned out 
to be true, it would suggest that in the gluon-mass description of QCD the ghost sector 
may not play such a central role as in the "ghost-dominance" picture 48|, obtained when 
working in the (conventional) Landau gauge. A preliminary study of these issues is already 
underway, and we hope to report its results in the near future. 
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APPENDIX A: MODIFIED ANGULAR APPROXIMATION 



In this Appendix we discuss the technical details related to the angular approximation 
and its modification. 

It is convenient to introduce the following quantities, appearing on the r.h.s. of (13.31) : 
I (q 2 ) = q 2 J d*kA (k)A(k + q), 
h(q 2 ) = q 2 J d 4 kA{k)A{k + q), 
I 2 (q 2 ) = J d A kk 2 A(k)A(k + q), 

I 3 (q 2 ) = iy d A kk 2 [k 2 -(k + q) 2 ]A(k)A(k + q). 



(A-l) 



Note that 



d A k 



{k-qf 



1 



1 



A(k)A(k + q) = -I 1 (q 2 ) + -I 3 (q 



(A-2) 



q" 4 K ' 2 

The angular integration of Io(q) may be carried out exactly, after shifting the integration 
variable k + q — > k, 



1 1 r fk A{k) 



7l~ 



(k + q) 2 
dyyA(y) + x I dyA(y) 



(A-3) 



For the remaining integrals we will use the angular approximation, appropriately modified 
to account correctly for their contributions to the running of the mass m 2 (x), i.e. terms that 
will enter in Eq. ( 14. 8ft . 

The standard angular approximation amounts to 



IT 



e(x-y)f(x) + B{y-x)f(y) 



(A-4) 



with z = x + y + 2 v /xy cosx, and 9(x) is the Heaviside step function. 
Defining 



(A-5) 



and using the superscript "A" to indicate that the aforementioned standard angular approx- 
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imation has been employed, we obtain 

PX poo 

If{x) = xA(x) dyyA(y) + x / dyyA 2 (y) , 

JO Jx 

I^x) = A(x) f X dyy 2 A{y)- fdyy'A^y) , 
Jo Jo 

Ife) = ^1 f X dyy\y- x )A{y). (A-6) 
x Jo 

Let us now check the faithfulness of these expressions for a special form of the massive 
propagator. Specifically, substitute in (1A-6I) A by the A of ( 14. 21) . apply the identity 
yA(y) = 1 — m 2 (y) A(y) , keeping only the terms linear in m 2 (y) and set A(x) — > to 
obtain 



f oo 

If(x)\ m2 = - I dym 2 (y)A(y)-x I dy m 2 (y)A 2 (y) 



o 

x ~ 1 , 



l£(x)\ m 2 = 2 dym 2 (y)A(y)-- dy y m 2 (y) A(y) , 
Jo x Jo 

I*(x)\ m2 = - / dyym 2 (y)A(y) - - / dyy 2 m 2 (y)A(y) . (A-7) 
^ Jo x Jo 

Now the result in (1A-7I) is to be compared with the direct calculation of Ii(x), before 
resorting to fl A-4I) . We begin with hiq 2 )] by substituting A by A into hiq 2 ) of (lA-lj) and 
isolating again the mass terms, we have 

I 2 (q 2 ) = Jd 4 kA(k) - J d 4 km 2 (k)A(k)A(k + q) , (A-8) 

and 

I 2 {q 2 ) - J 2 (0) = J d 4 km 2 {k)A{k) ^A(k) - A(k + q)j . (A-9) 
This is the exact mass dependence. To this last expression we now apply ( 1A-4I) . to get 

lf{x)\ A = f X dyym 2 (y)A(y)\A(y)-A(x) 
Jo L 



X 1 PX 



dym 2 {y)A{y) - - I dy y m 2 (y) A(y) + 0{m 4 ) . (A-10) 
'o x Jo 

Comparing 7™ (x)\ A with 7 2 A (»| m 2 [Eqs.( fA^T0|) and (fATfl) we see a discrepancy of a factor of 
2 in their first term. The simplest term, quadratic in A, that could correct this discrepancy 
is | JL X dyy 2 A 2 (y); adding it to the result of the standard angular approximation, I£(x) 
[Eq. flA-61) ]. will lead to a modified 7 2 4 (x)| m 2, that will coincide with I™ 2 (x)\ A . 
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We will next repeat the same exercise for h(q 2 ). It is elementary to show, by substituting 



A -> A into h(q 2 ) that 



I 3 (q 2 ) = 1 J d 4 k [m 2 {k + q)- m 2 {k)] k 2 A{k)A{k + q) + J d 4 kA(k) , (A-ll) 
and therefore, applying the angular approximation and subtracting at x — 0, 

If(x)\ A = ^ [ X dyy 2 [m 2 (x)-m 2 (y)} A(y) 



x 

2/„\ 1 /-a; 

2™2 

1 till II 



m (x) 1 
~~ 2 x 



dyy 2 m 2 {y)A{y). (A-12) 



o 



This term is subleading for either type of running of the mass (logarithmic or power- law). 
This important property is obviously not captured by the approximate expression of I£(x)\ m 2 
in (1A-7I) . since the first term would furnish (erroneously) a leading order contribution in 
the case of power-law running. This shortcoming may be remedied by simply adding to 
J3 (x), Eq. (lA-6p . the term i f* dyy 3 A 2 (y); thus, the subleading nature of the term J 3 (g) is 
preserved. 

We finally turn to Ii{q 2 )- The most immediate way to see that the standard angular 
approximation for I\, given by the first equation in (1A-6[) . mistreats the running of the 



mass, is to set A(y) — > A(y) and then m 2 (y) = m 2 in the initial expression for I\{q 2 ) 
(second equation in (1A-1I) ). Then h(q 2 ) gets reduced to a standard one- loop integral, 

/d 4 k 
( p + m2) p + g)2 + m2] ' ( A " 13 ) 



and thus h(q ) is given by 



W) = -q 2 



, , m 2 \ ( D + 1 

1/2 



(A-14) 



where c — — f + 7 — ln47r and D = 1 1 + J . For large values of q 2 we have that the 
finite part of h(q 2 ) goes like 

Ii(x) ~ — x lnx — 2m 2 lnx . (A-15) 

Notice that the last term, which determines the contribution to the running of the mass, 
comes with weight 2 compared to the logarithmic term determining the asymptotic running 
of A _1 (x) (remember that I\{q 2 ) is to be multiplied by (—6), so the first term of (1A-I4j) 
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provides the RG logarithm). Instead, setting A(y) — > A(y) and m 2 (y) = m 2 in J-j 4 (x) and 
Ii(x)\ m a, [Eqs.( lA-6l) and ( 1A-7I) ]. we find as leading contribution 

Ii(x)\ m 2 ~ — m 2 lnx 

Ii(x) ~ — xlnx — m 2 \nx. (A-16) 

Unlike the (exact) expression of flA-15j) . now the two logarithms appear with equal relative 
weight. Evidently, the angular approximation captures correctly the leading (RG) logarithm, 
but furnishes only half of the m 2 In x contribution. 

There is another way to verify that the correct relative weight between the two logarithms 
of h{q 2 ) is 2. When the spectral representation for A(q) is used, 

A(? 2 )= jdX 2 / ( . A 2 V , (A-17) 
J q z — X z + ie 

the "kinetic" term of the resulting (linearized) SDE, i.e. the exact analogue to —Ii(q 2 ), is 
given by 

B(q 2 ) =q 2 Jd\ 2 p(X 2 ) j d 4 k A(jfc) A(k + q) , (A-18) 
and may be easily cast in the form 

B(x) =cx + xj^ 1 dy (l-^j A(y), (A-19) 

where c is a (divergent) constant, to be absorbed into the wave-function renormalization of 
the gluon field. When determining the contribution of B(x) to the running of m 2 for large 
x, it would be wrong to simply set (l — -j) 1 ^ 2 — ► 1; instead, one must expand this term to 
first order. Setting A(y) — > A(y) we obtain 

px/i px/4 

B(x) = cx + x / dyA(y) — 2 / dyyA(y) 
Jo Jo 

nx/4 

~ c'x + xlnx + 2 / dym 2 (y)A(y) , (A-20) 
Jo 

where d includes now some irrelevant (finite) constant term. Again, the mass term is 
multiplied by a factor of 2 compared to the RG logarithm. 

The simplest way of remedy this discrepancy is to modify the result of the angular 
approximation for h(q 2 ), adding the term | J* dyy 2 A 2 (y); note that this is exactly the term 
that has been added to hiq 2 )- 
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Thus, finally, we arrive at the following modified expressions, which capture correctly the 
leading m 2 -dependence of the SDE: 
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